First limit on inclusive B ^ XgUu decay 
and constraints on new physics 

(Erratum) 

Yuval Grossman", Zoltan Ligeti^ and Enrico Nardi'^ 

""Department of Particle Physics 
\ Weizmann Institute of Science, Rehovot 76100, Israel 

California Institute of Technology, Pasadena, CA 91125 

(July 16, 1996) 



=3 

1—5 



O 

in 

ix. 



X 



In a recent article we presented the derivation of the first hmit on the inclusive 
B ^ XgV v decay rate. Our work stemmed from the observation that a published ALEPH 



> 

' bound on BR(i? r u) 0, inferred from the absence of large missing energy events in B 



decays, implies also a limit on BR(i? ^ XgUu). We estimated this limit by comparing the 
missing energy spectra in the two decay modes. Theoretically, i? — z/ z/ is a very clean 
process, which is also sensitive to several possible sources of new physics |T|. Therefore, a 
dedicated experimental search is important. The result of such an analysis by the ALEPH 
Collaboration, using the full LEP-I data sample, is to appear soon [§. 

As a result of discussions with members of the ALEPH Collaboration [^, we found two 
errors in our numerical results. The corrected analysis yields a limit weaker than our original 
result [in by about a factor of three. The corrections included here are formally beyond the 
free quark decay result at tree level (of order a^, Aqcb/ rrih, and higher). However, they 
affect the final limit significantly, due to the specific details of the experimental analysis. 

First, due to a mistake in our Monte Carlo code, the energy of the Xg hadronic final 
state in the B rest frame was evaluated as a fraction of the b quark mass, rather than as a 
fraction of the B meson mass. This resulted in a harder missing energy spectrum than the 
true one. Correcting for this weakens the limit by about 35%. 

Second, we neglected the effects of the invariant mass distribution of the final hadrons, 
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which we (imphcitly) assumed not to extend to high mass states. However, the invariant 
mass of the Xg system can be well above 1 GeV Because of the large boost into the 
LEP laboratory frame and the high missing energy range (> 35 GeV) used in the analysis 
(which is large compared to the average B meson energy at LEP), neglecting this effect 
yielded a missing energy spectrum in the laboratory frame considerably harder than the 
correct one. 

It is not straightforward to include properly this second effect into the analysis. However, 
an estimate can be obtained by approximating the invariant mass spectrum of the Xs system 
with a Gaussian distribution with mean (yu) and variance (a) fitted to the averages (M^) 
and (M^), as given in 0. Such a fit yields /i ~ 1.35 GeV and a ~ 0.6 GeV. There is about 
a 10% uncertainty in these values, due to their dependences on the heavy quark effective 
theory parameter A (~ itlb — rrih) . Treating the Xg mass distribution as independent of the 
energy spectrum, and including a theoretical uncertainty of about 15% related to the fitted 
values of fi and a, we find that the original limit given in is weakened by about a factor 
of three. While small values of Ex favor small values of Mx (beyond the trivial kinematic 
constraint Ex > Mx), the correlation between Ex and Mx only slightly improve the bound. 

In conclusion, the bound on BR(i? ^ XgUu) given in the Abstract and in Eqs. (1.2) 
and (6.1) of [IH is weaker by about a factor of three, while the limits on the new physics 
parameters, collected in Table 2, are weaker by about a factor 1.8. Until the ALEPH bound 
on the B ^ XgU u branching ratio will appear, the limit BR(i? XsUu) < 1.3 x 10^^ 
should be used, instead of that quoted in [|l|. 
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Abstract 

The inclusive B Xs v v decay rate, on which no experimental bound exists 
to date, can be constrained by searching for large missing energy events in 
B decays. Carefully examining the experimental and theoretical aspects of 
such an analysis, we argue that the published ALEPH limit on BR(i? tv) 
implies, conservatively, the bound BR(i? Xg ui?) < 3.9 x 10""^, which is less 
than one order of magnitude above the standard model prediction. The LEP 
collaborations could significantly improve this bound by a dedicated experi- 
mental analysis. We study the constraints this new limit imposes on various 
extensions of the standard model. We derive new bounds on the couplings of 
third generation fermions in models with leptoquarks, and in supersymmetric 
models without R-parity. We also constrain models where new gauge bosons 
are coupled dominantly to the third generation, such as TopColor models and 
models based on horizontal gauge symmetries. For models which predict an 
enhanced effective bsZ vertex, the constraint from B Xg is competitive 
with the limits from inclusive and exclusive i? — > decays. 
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I. INTRODUCTION 



Recent progress in experiment and theory has made flavor changing neutral current 
(FCNC) B decays a stringent test of the Standard Model (SM) and a powerful probe of 
New Physics (NP). The CLEO Collaboration observed the exclusive decay B —>■ K* 
as well as the inclusive decay B Xs'~f The UAl upper limit on the inclusive decay 
B ^ Xg ji^ ji^ 0, and the recent CLEO |^ and CDF upper limits on the exclusive 
decays B — > are less than one order of magnitude above the SM predictions. 

These decays are likely to be observed within the next few years. These experimental 
results constitute a set of very strong constraints on several possible sources of NP. 

The FCNC decay B XgUv is also very sensitive to extensions of the SM, and provides 
a unique source of constraints on some NP scenarios which predict a large enhancement of 
this decay mode. In particular, the B ^ XgVj. mode is very sensitive to the relatively 
unexplored couplings of third generation fermions. However, no experimental upper bound 
on this decay mode has been established to date. 

The decay B XgVv can be searched for through the large missing energy associated 
with the two neutrinos. Using such techniques, the ALEPH L3 |p, and OPAL 

collaborations have been able to measure the inclusive B ^ X^-tv decay rate [|10|. The 
large missing energy in this case is associated with the two neutrinos in the decay chain 
B ^ XcT u followed hj r ^ u X. From the absence of excess events with very large missing 
energy, ALEPH also established the 90%CL bound 

BR(B ^ rz/) < L8 X 10-^ (1.1) 

In this paper we point out that a similar analysis of the same data also implies a limit 
on BR.{B ^ XgUu). While a detailed and complete analysis of this set of data can only 
be performed by the ALEPH Collaboration, we show that using some conservative and 
simplifying assumptions it is possible to derive the bound 

BR(5 ^ X,z/z/) < 3.9 X 10~\ (1.2) 
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which is less than one order of magnitude above the SM prediction |TT 



BR^"(5 X,z/z/) ^ 5 X 10 
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(1.3) 



We expect that a dedicated analysis of the LEP collaborations will strengthen our bound. 
If background subtraction can be performed, or extra experimental cuts can reduce the 
background, then we estimate that using the full LEP-I data sample a 90%CL bound of 
order 



could be within the reach of the LEP experiments. 

In section II we discuss, in a model independent manner, the inclusive B ^ X^v v decay 
rate and the missing energy spectrum. In section III we describe the theoretical issues 
involved in relating the limits on large missing energy events to the B ^ XsV v decay rate, 
and we estimate the corresponding theoretical uncertainties. This section also contains the 
details of the derivation of our bound ( |1.2D . While the limit we obtain does not allow a 
direct test of the SM, it still implies stringent constraints on various new physics scenarios. 
In section IV we analyze the implications of existing measurements of FCNC processes, 
particularly B ^ Xs'y and i? — yU+ yU^, for 5 ^ z/ z/, in various classes of NP models. 
Section V contains an extensive discussion of the constraints implied by the bound ( |1.2|) on 
various extensions of the SM. We derive new constraints on models in which the couplings 
of the third family fermions differ from those of the first two generations, and we present 
numerical limits on several NP parameters. Finally, section VI contains a summary of our 
results and the conclusions. 

II. THE B ^ XsvD DECAY RATE 

Prom the theoretical point of view, the decay B XsV v is a, very clean process, since 
both the perturbative and the non-perturbative 1/m^ corrections are known to be small. 



BR(S ^ X, z/ z/) < (1 - 2) X 10 
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(1.4) 
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Furthermore, in contrast to the decay B ^ Xgi^ i~, which suffers from background such 
a.s B —>■ X s J / ip Xgi"^ i~, there are no analogous long-distance QCD contributions, since 
there are no intermediate states that can decay into a neutrino pair. Therefore, the decay 
B ^ Xgiy iy is well suited to search for and constrain NP effects. 

As our aim is to derive constraints on NP scenarios, we discuss the missing energy spec- 
trum in a model independent framework. Limits on NP parameters can then be derived by 
comparing the experimental bound with the theoretical expressions as derived in specific 
models. A model independent expression for the missing energy spectrum can be straight- 
forwardly obtained from the general result for muon decay [12| . Under the only assumption 



of two component left-handed neutrinos (possible neutrino mass effects are at most of or- 
der iriy/mi, < 10^^), the most general form of the four-fermion interaction responsible for 
B — s> Xq Vi Vj reads 

C = ClOl + CrOr, (2.1) 

where 

Ol = [qL 1, &l] K YK] , On = [qn b^] [vl 7^^i] ■ (2-2) 

Here L and R denote left- and right-handed components, q = d,s, and i,j = e,fi,T. In 
this article we adopt the notation that a generic B meson contains a b quark, rather than 
a b quark. As the flavors of the decay products are not detected, in certain models more 
than one flnal state can contribute to the observed decay rate. Then, in principle, both Cl 
and Cr carry three indices q, i, j, which label the quark and neutrino flavors in the flnal 
state. Throughout our discussion, we shall only keep track of these indices when they are 
important, otherwise we will suppress them. 

At lowest order, the missing energy spectrum in the B rest-frame is given by 



l£(^^|iilM = i!l(|C,P + |C„P)5('-.-)- (2-3) 

Here we have not yet summed over the neutrino flavors. The function S{r,x) describes the 
shape of the missing energy spectrum 



S(r, x) 



[1 — x) (4x — 1) + r (1 — 3a;) — Gr/v^ (1 — 2x — r) 



(2.4) 



The dimensionless variable x = Emiss/^b can range between (1 — r)/2<x<l — ^Jr, and 
r = ml/ml. The parameter 77 = — Re(CL C]^)/(|Cj^P + |C/jp), which is the analog of the 
Michel parameter in //-decays, ranges between —\<rj<]^. 

In the SM, B ^ XgU u proceeds via W box and Z penguin diagrams, therefore only Ol 



is present and r] = 0. The corresponding coefficient reads |T3[] 

V2GFa 



^SM 



71 sin^ 9]y 



V,lVtsXo{xt) 



(2.5) 



where Xt = rn^/m^, and 



X.(x) 



X 



2 + X 3x — 6 

+ 7 7T7T InX 



X 



(2.6) 



1 ■ (X-1)2 

In the limit of large top quark mass, Xq has a quadratic dependence on rrit, Xo{xt) ~ Xt/8. 
Therefore, the main source of uncertainty in the SM prediction for the total decay rate 
comes from the uncertainty in m^. 

The leading 1 / ml and corrections to the SM result (calculated in the free quark decay 
model) are known. The as correction to the total decay rate is given by replacing Xo{xt) in 
Eq. (O) by 



Xoixt) 



a. 



Xoixt) + X,{xt) 
47r 



3tt\ 4 / 



(2.7) 



Here the second term represents the correction to the matrix element of Ol- This term can- 
cels almost completely when the B ^ X^v v branching ratio is normalized to the semilep- 
tonic B ^ XcCv rate (see Eq. (|2.9|) below). The first term contains the QCD correction to 
the box and penguin diagrams. We do not display here the explicit form of Xi(xt), which 



can be found in |11|. The most important effect of this correction is to reduce the scale 



dependence of the SM prediction from about ±10% to below ±2% |]TT|. 

The correction to the contribution of the Ol operator to the missing energy spec- 



trum can be read off from The result is the following modification of the function 

S{r,x) in Eq. (U) 



iS(r, x) S{r, x) 



^ [2(1 - x? (19 - 38x + lOx^) - r(61 - 23r - 122x + 52x^)1 

. ^,^2 _ ^ IQmi 



1 — xY — r 

[(1 - x)(47 - 126x + 96x2 _ 20a;3) _ ^(yg - 23r - 125x + 52x^)1 

2m 



+ 



(l-r)3(5(i(l-r) -X 



l-r)^(5' -(l-r)-x . (2.8) 



Here Ai and A2 are related to the kinetic energy of the h quark inside the B meson and to 
the mass sphtting between the B and the B* mesons, respectively. Experimentally, A2 ~ 
0.12 GeV^, and following the discussion in [|1^, we use < — Ai < 0.5 GeV^. When integrated 
over the spectrum, this correction amounts to about 3% suppression of the total decay rate. 
Even this small correction cancels almost entirely when the B ^ XgV v branching ratio is 
normalized to the semileptonic B ^ XcCv rate (see Eq. ( |2.9| ) below). As we will discuss 
in the next section, although the B —>■ Xg u u branching fraction in the SM is known rather 
precisely, the theoretical prediction for the missing energy spectrum is more uncertain. 
In several NP models Or is also present. The structure of the operator product ex- 



pansion [|T6|-|T^ shows that to all orders in the and l/m^ expansions, both the per- 
turbative and non-perturbative corrections to the contribution of On to the missing en- 
ergy spectrum are identical to those of Ol- This holds as long as the phase space is 
symmetric in the two leptons (this can be violated only by negligible neutrino mass ef- 
fects). Thus, the shape of the missing energy spectrum is model independent, up to possible 
small interference effects between Ol and Oji, of order ms/rrih. Once Cl and Cr are com- 
puted in any particular model, the unknown contributions to the total decay rate are only 
0{al; asms/nib] A^/m^; A^/m^; A^/m^), where A denotes some scale of order Aqcd- 
In some NP models Ol and Opt are simultaneously present with coefficients of comparable 
size, giving rise to interference effects proportional to 1], which modify the shape of the 
missing energy spectrum. Due to the ms/mf, suppression, this effect is always small, except 
close to the endpoint region (x ~ 1 — y/r), where the leading term in ( |2.4|) also becomes of 



order y^. Therefore, close to the endpoint, the corrections to the shape of the spectrum 
could be relevant. Thus, besides the quark mass ratio, an additional uncertainty in ( p.3| ) 
is associated with the value of rj. The softest missing energy spectrum, and thus the most 
conservative bound on the B ^ XsV v branching ratio, is obtained by using a large quark 
mass ratio r ~ 0.002 (corresponding to ^ 0.2 GeV and m;, ^ 4.8 GeV), and r] = The 
bound on the branching fraction derived using these values of r and rj holds in any model. 

In order to find the constraints on NP, we will express the total decay rate in terms of two 
"effective" coefficients C l and Cr, which can be computed in terms of the parameters of any 
NP model and are directly related to the experimental measurement (see (|2.12| )). To remove 
the large uncertainty in the total decay rate associated with the ml factor, it is convenient to 
normalize BR(i? — > z/ z/) to the semileptonic rate BR(i? XcCu), since the experimental 
value of the latter is known quite precisely. The contribution from B X„ e z/, as well as 
possible NP effects on the semileptonic decay rate are negligible. In constraining NP, we 
can also set = and neglect both order and corrections. This is justified, since 

when averaged over the spectrum these effects are very small, and would affect the numerical 
bounds on the NP parameters only in a negligible way. For the total B Xg z/j i?j decay 
rate into all possible q = d,s and i, j = e , fi ,t final state fiavors, we then obtain 

BRjB^Xui.) ^ Cl + Cl 

BR{B ^X^ev) fps{ml/ml) ' ^ ' ^ 

where fps{x) = 1 — 8a; + 8x'^ — x^ — 12x^ In x is the usual phase-space factor, and we defined 



cl = ^T.\cr\ c'r = ^ ^- (2.10) 



The SM prediction for the branching ratio in Eq. ( |1.3| ) is obtained by inserting into 
Eq. (U) the semileptonic rate BR{B ^ X^eij) 10.5% fps{ml/ml) 0.5, together 
with 

(C'^)'^3.8xlO-^ (2.11) 

which follows from (^ by using \V;*^Vts\/\Vcb\ ^ 1, \Vcb\ ~ 0.04, a{mz) ^ 1/129, sin^ ^vk ~ 
0.23, and rrit ^ 180 GeV [|l^]. In a SM analysis it would be more natural to factor out 



from the definitions of Cl.r the small mixing angles V^lVts as well as a, resulting in a 
dimensionless coefficient of order unity. However, there is no reason to do so for Cr, and 
since in several NP models C^jj are induced at the tree level, they do not contain any small 
parameters analogous to a or to the CKM angles. In fact, often the only suppression factors 
in these coefficients come from inverse powers of some large mass scale. 

Comparing the limit ( |1.2| ) with ( p.9| ) yields the following constraint on possible NP 
contributions: 

Cl + Cl< 3.0 X 10"^ 

In section V we will constrain NP models by comparing this limit with the various theoretical 
predictions for the coefficients Cl,/j. We will quote the bounds that can be derived on each 
NP parameter, even in those cases when Ol and Or are simultaneously present and certain 
combinations of the NP parameters may be better constrained. 

III. EXPERIMENTAL ANALYSIS 

In this section we discuss the theoretical uncertainties involved in the experimental anal- 
ysis at LEP, and we derive the bound ( |1.2| ). To obtain a limit on B ^ XsV v using the data 
from ALEPH , we have to estimate the relative efficiency oi B X^v v and B ^ rv io 
pass the experimental missing energy cut (ALEPH used the cut -Emiss > 35GeV). We use 
the B ^ XsV v missing energy spectrum as given in Eqs. ( p.3|) and (p.4|) , together with the 
conservative values r = 0.002 and r/ = — | which, as discussed in the previous section, lead 
to the softest missing energy spectrum. As a result, the bound we derive holds in any model. 

A. Theoretical uncertainties 

The largest theoretical uncertainty is related to the reliability of the theoretical calcu- 
lation of the missing energy spectrum near its endpoint (a; ~ 1 — \/r). While both the 



BR(5 ^ Xsuu) 
3.9 X 10-4 



(2.12) 
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and corrections to the total decay rate are small, this is not true point-by-point for 

the differential energy spectrum. 

The l/nib expansion for inclusive semileptonic decays of hadrons containing a single h 
quark based on an operator product expansion |T6[ is not reliable when the invariant mass 
of the decay product hadronic system is small. In the endpoint region, near x = 1 — i/r, 
there are large corrections to the free quark decay prediction, and the spectrum has to be 
smeared to get a reliable result. This smearing region has to be chosen large enough so 
that after smearing, the corrections in ( p.8|) produce only small effects on the decay 



spectrum [jT^ • Following this criterion, we find that in the present case the smearing region 



extends about 0.5 — 0.7 GeV from the maximal value of the missing energy. 



The order corrections to the missing energy spectrum can be read off from pO[, and 
are less problematic. While there are large logarithms for small values of the hadronic 
invariant mass, the perturbative correction to the hadron energy spectrum is smooth in the 
limit of vanishing hadronic energy. Thus, the perturbative corrections to the missing energy 



spectrum are small and reliably calculable in the region relevant for our analysis pO| , P |. 

The large non-perturbative corrections near maximal missing energy could represent 
a problem for the analysis of B ^ XgU u, since the experimental cuts select events with 
large missing energy, thus enhancing the weight of the problematic region. However, the 
experimental search is at the endpoint region in the laboratory frame. Events with large 
missing energy in the laboratory frame are not necessarily those with missing energy close 
to the endpoint in the B rest-frame. We need to take into account the large boost from the 
B rest-frame to the laboratory frame. We define Pe[GeV](-£'cut[GeV]) as the fraction of the 
events with -Emiss > Ecut in the laboratory frame, that come from the missing energy region 
in the B rest-frame between the endpoint and Eq = mi, — rris — e. To simulate the B hadron 
energy in the laboratory frame, we use in our Monte Carlo code the Peterson fragmentation 
function and nib ~ 4.8 GeV. We find the following representative numbers 
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Po.5(40) ^0.12, Po.5(35) ^0.11, Po.5(30) ^0.09, 

Po.7(40) ^ 0.21 , Po.7(35) ^ 0.19, Po.7(30) ^ 0.17 . (3.1) 



The results in (|3.1|) show some sensitivity to the size of the smearing region, while they 



suggest that varying the experimental cut has a much smaller effect. We have also checked 
that these figures are not very sensitive to the choice of the quark masses and of 1]. 

We conclude that the missing energy spectrum in the laboratory frame of the LEP 
experiments can be estimated reliably from the spectator model. However, the uncertainty 
related to the endpoint region of the missing energy spectrum in the B rest-frame should 
be taken into account. To estimate this uncertainty, we rely on the fact that the size of the 
endpoint region is chosen such that the integration over it can be trusted. Assigning to the 
fraction of the events coming from endpoint region the lowest possible value of the missing 
energy, — — e, we can bound the uncertainty related to using the spectator model. 
For e = 0.7 GeV we find in this way a 5% reduction in the total number of events that pass 
the 35 GeV cut. We consider this as a reliable estimate of the theoretical uncertainty of the 
analysis related to the use of the spectator model. 

A final remark is in order. At LEP, b quarks hadronize into a variety of b hadrons. 
However, by the time it decays, the b quark is contained either in a P meson or in a A^ 
baryon. The fraction of A;, baryons at LEP has been measured to be about 10% |2^. In the 
limit when the b quark is treated as very heavy, the missing energy spectrum from B decays 
should be similar to that from Af, decays, up to small effects originating from the polarization 
of the baryons. A possible indication for significant heavy quark symmetry violating effects, 
is the experimentally measured lifetime ratio between the B and the Af, hadrons, which 
appears to be larger than the theoretical prediction (for a recent discussion see, e.g., [p3[] ). 
However, the resolution of this problem is most likely related either to experimental issues 
or to the theoretical calculation of the hadronic decay widths, so it is probably irrelevant 
for the analysis in this paper. The polarization of the Af, baryons produced at LEP has 
been measured to be —30 ± 30% [P5|. A non- vanishing left-handed polarization enhances 
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the missing energy from Af, decays, and therefore neglecting it is conservative. 

B. Derivation of the bound 

After this discussion, we are ready to estimate the bound on BIi{B —>■ Xg u u). To 
translate the ALEPH bound on B ^ t u ^ into a limit on B ^ XsV we need to compare 
the theoretical predictions for the fraction of events that pass the missing energy cut for the 
two decay modes. The resulting bound on B ^ XgUv is stronger than that on B ^ tv for 
the following reasons: 

{i) The B TV decay is allowed only for the charged B mesons, while all h flavored 

hadrons can decay through the parton level process b —>■ su i?; 
(ii) In order to reject background from semileptonic B decays, only the hadronic r decays 

were used in the ALEPH search for B t z/.[] 
{in) The missing energy spectrum is somewhat harder in B ^ XsV v than in B ^ t 9 

decays, increasing the efficiency of the analysis. 
To evaluate the last factor [iii), we need to estimate the missing energy spectrum for the 
B ^ TV decay as well. In estimating the t ^ v X missing energy spectrum, we consider 
only two-body (20%) and three-body (80%) hadronic r decays. Since a significant fraction 
of hadronic tau decays are four- and five-body decays, the resulting missing energy spectrum 
for 5 ^ r z/ is harder than the real one. Therefore we get a conservative upper bound on 
the B XgV i? branching ratio, i.e., weaker than what a detailed analysis would obtain. 
We take into account that the r polarization decreases the efficiency of the B t u analysis 
by about 20%, as estimated in the ALEPH analysis 0. We find that in our approximation 
the efficiency of the B —>■ XsV v decay to pass the ALEPH cut -Emiss > 35 GeV is about 15% 
larger than that oi B ^ tv followed by hadronic r decay. 



*We thank Ian Tomalin for pointing this out to us. 
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Collecting the results of our previous discussion, we can finally estimate that the bound 
on BR(i? —>■ XgU 1?) is stronger then the ALEPH bound on BR(i? ^ t u)hj an overall factor 

R = Rb± i?hadr ReS ~ 0.21 . (3.2) 

Here Rb± ~ 0.37 is the ratio of mesons to the total number of b hadrons at ALEPH 



2^ and -Rhadr ~ 0.65 is the hadronic r branching fraction [^. The factor R^s accounts 
for the efficiency of the missing energy cut in the B ^ tv decay followed by hadronic r 
decays, relative to that in B ^ XsV v decays. We conservatively use the upper bound 
i?efr = 0.90 which includes the 5% theoretical uncertainty related to the reliability of the 
spectator model spectrum. We expect the uncertainty related to the use of the Peterson 
fragmentation function to cancel to a large extent from the estimate of the relative efficiency 
of the B ^ XsV v and B ^ tv analyses. Using ( p.2| ) and the experimental bound ( |1 . 1|) , we 
find the limit given in Eq. ( p..2|) . 

We would like to emphasize that the derivation of our bound relies on a set of conservative 
simplifying assumptions, which could be avoided in a dedicated experimental analysis. Our 
limit ( |1.2| ) should be considered as a conservative upper bound, to which we purposely do 
not (and cannot) assign a confidence level. We hope that a more complete and detailed 
investigation by the LEP collaborations will be carried out. The reward of such an analysis 
could be a bound of order (1 — 2) x 10"^, that is only a factor 2 — 4 above the SM prediction. 



IV. CONSTRAINTS B ^ Xsvv FROM OTHER PROCESSES 

Models that can give rise to large new contributions to the B ^ XgV v decay often 
predict an enhancement of other FCNC processes as well. In this section we analyze what 
constraints the existing experimental data on FCNC processes imply for various NP models. 
For each specific model, these constraints result in upper limits on the allowed B Xgiy u 
decay rate. Bounds on BR.{B ^ XsV v) can be derived from the limits on rare processes, 
such as — s> K — > vri/z/, e^^, K — K and B — B mixing. The most restrictive 
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constraints are imposed by the measurements of the radiative decay B ^ Xg'j, and by the 
hmits on inclusive and exclusive B Xgi^ i" decays. 

The radiative B —>■ Xs'~f decay proceeds via photon penguin diagrams, and therefore it is 
not directly related to i? ^ z/. However, in many models the details of the underlying 
physics imply relations between the Z and the photon penguins. In all these models the 
recent CLEO measurement 

BR{B ^ X, 7) = (2.32 ± 0.51 ± 0.32 ± 0.20) x 10"^ (4.1) 

which is in agreement with the SM, forbids large deviations from the SM prediction for the 
B ^ XsV v decay rate as well. 

On the other hand, a large class of NP models predict (or can accommodate) an enhanced 
hsZ effective vertex without giving rise to a large enhancement of the 657 effective coupling. 
Then the constraints from inclusive and exclusive B Xsi~^ i~ decays are important, as 
these decays, like B ^ XsV v., are dominated by Z exchange. In these models, a naive 
estimate of the ratio of inclusive rates gives BR(i? — > z/ z/)/BR(i? Xgi^ i^) ~ 6. The 
factor of six enhancement arises due to a factor of approximately two in the ratio between 
the neutrino and the charged lepton couplings to the Z, and a factor of three from the sum 
over the neutrino flavors. A more precise calculation which includes the photon exchange 
contribution to B ^ Xs i~ bounded by the CLEO measurement of i? ^ X^ 7, as well as 
the sizable QCD corrections to B ^ Xsi^ i~ p5|j26[| , can increase the above ratio up to 7. 
Hence, for this class of models, the UAl experimental limit on the inclusive B ^ Xg fi^ fi~ 
decay iQ 

''^The UAl experiment searched for events in the region 3.9 < E"^^ < 4.4 GeV. However, the theo- 
retical prediction for the spectrum is uncertain in this endpoint region. Moreover, the theoretical 
spectrum shown in Ref. |^ (and presumably used to relate the bound on the number of events 
in the experimental window to the quoted bound on the total branching fraction) seems to be in 
disagreement with that calculated in [|^. Despite these uncertainties we use the published bound. 
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BR(5 ->X,/i+/i-) < 5 X 10-^ (4.2) 

implies 

BR{B -^X,vv)< 3.5 X 10"^ =^ c\ + c\< 2.7 x 10"^ . (4.3) 



This limit is comparable with the limit ( |1.2| ), however it is weaker than the expected LEP 
sensitivity ( |1.4| ). This underlines the importance of a more detailed experimental analysis 
aimed at searching for B ^ XgV and shows that LEP measurements could compete with 
future new data from CLEO and CDF in constraining models which predict an enhanced hsZ 
coupling. Moreover, as we have emphasized, the neutrino mode is particularly interesting 
since it is theoretically cleaner than the charged lepton modes: (z) there are no long-distance 
QCD effects; {ii) the short-distance QCD corrections are small; {Hi) there is no photon 
penguin contribution, and therefore this process can be straightforwardly related to the 
effective hsZ vertex. In conclusion, although the constraints resulting from our limit ( |1.2|) 
are numerically slightly weaker than those implied by S ^ Xsfi~^ fi~ ( [4.3|) , we think that 
they are more reliable. 

The exclusive dilepton decays are also sensitive to an enhanced bsZ vertex, and provide 
additional constraints. The exclusive decay modes B K^*^ i~ have been searched for 
by the CLEO |^ and the CDF ^ collaborations. For our purposes, the most restrictive 
limit has been established by CLEO 

BR{B ^ K* e+ e~) < l.Q X IQ-^ . (4.4) 

However, the interpretation of this limit is obscured by the significant model dependence of 
the exclusive B K^*^ form factors. The estimates for the ratio 

BR(5^K*e+e-) ,^ , 

^" BR(i^^X.e^e-) ' ^'''^ 



range between 0.10 and 0.35 |2^. Clearly, the value of p is crucial for relating the limits 



on exclusive decays to a limit on the FCNC transition at the quark level, for which the SM 
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predicts [p5[] f1 

BR(5 ^X,e+e-) ^ 10 X 10-^ BR(S ^ X, /i+ /i") ^ 7 x 10-^ (4.6) 

It is also questionable whether the above estimates of p, based on SM computations, are 
reliable for the analysis of NP scenarios. In fact, in the SM a significant contribution 
to -B ^ K* e"*" e~ comes from soft photons, corresponding to small e^e" invariant mass. 
However, in the NP models under consideration Z penguins dominate, and therefore p is 
expected to be smaller than in the SM. 

The only known model independent way to predict the B K^*^ form factors (besides 
lattice calculations) is to relate them to measurable semileptonic decay form factors using 



heavy quark symmetries ||28|. However, until experimental information on the B ^ piu 
decay spectrum becomes available, the dilepton spectrum can only be predicted in the small 
window 4.0 < rrie+e- < 4.4 GeV from the D —>■ K* i'^u data. (Even after the measurement 
of the B —>■ piu spectrum, the unknown symmetry breaking corrections will amount to an 
uncertainty of order 20 — 30%.) 

Rare K decays also imply bounds on Z penguins. However, in relating them to the 
B ^ XsV v decay rate large uncertainties arise from poorly known quark mixing angles. 
The only measured rate is Kl — > p^ p^ , which receives large long-distance corrections. 



Therefore, the short-distance parameters can only be extracted with large uncertainties P9 |. 
The decay ^ tt^ u i?, however, receives negligible long-distance contributions. While the 
existing experimental bound is about fifty times the SM prediction, an order of magnitude 
improvement is expected from Brookhaven in the coming years. In summary, at present, for 
all the models in our analysis, rare K decays are less constraining than the rare B decays 
discussed above. 

Once these considerations are taken into account, it appears that the most reliable con- 



•f-The spread in the SM predictions in the literature is in part due to the different signs for the 
O7 — Og interference term in various papers. (We use here the notation of Ref. [^.) 
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straints on an anomalous effective bsZ vertex, apart from B XsV y, come from the limits 
on the inclusive B ^ fi^ fi~ decay rate. Numerically similar bounds (or possibly slightly 
better - depending on the adopted values of model-dependent hadronic form factors) are 
provided by the bounds on exclusive B K^*^ t~ decay rates. 



V. NEW PHYSICS 

Following the discussion in the previous section, it is useful to classify the NP mod- 
els which could enhance the SM prediction for B ^ X^v v into three classes: (A) Highly 
constrained models; (B) Weakly constrained models; (C) Unconstrained models. 

The first class (A) includes models in which the existing bounds on other FCNC processes 
(mainly 5 ^ X,, 7) imply that the rate for B XgVV cannot exceed the SM prediction by 
any factor larger than two. Thus, a bound owB ^ X^vv oi the order of the LEP sensitivity 
( |1.4[ ) does not imply any new constraint on the underlying NP. For completeness, we briefly 
explain below why the B X^vv branching fraction must be close to its SM value in the 
most popular models belonging to this class: the minimal supersymmetric standard model; 
multi Higgs doublet models; left-right symmetric models. 

To class (B) belong all models which predict (or allow for) a large effective hsZ (or hsZ'^ 
vertex. We call these models "weakly constrained" , as the limits on inclusive and exclusive 
B ^ Xgi'^ l~ already constrain them. The limits from ( |1.2|) will not represent any numerical 
improvement over the existing bounds on the parameters of these models. However, our new 
limits are more reliable, since the theoretical uncertainties involved m B ^ XsV v are very 



small. Moreover, if a bound of the order ( p..4| ) can be obtained by the LEP collaborations. 



then the constraints from B ^ XsV v will become markedly stronger than the present ones 
from B ^ Xgl^ l~ . A list of interesting NP models belonging to this class includes: models 
with additional Q = — | isosinglet quarks; models with large hsZ' vertex; models with 
an anomalous effective tcZ vertex; models with a heavy fourth generation; models with 
anomalous WWZ couplings; a class of extended technicolor models. These models will be 
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discussed in some detail below. 

In the unconstrained models of class (C), the couplings responsible for enhancing 
B ^ XsV v are to a large extent independent of those constrained by any other existing 
experimental bound. Therefore, even a, B ^ XsV v decay rate orders of magnitude above 
the SM prediction is still consistent with the existing constraints, and the new limit ( |1.2[ ) 
represents the most stringent bound on the corresponding NP parameters. Examples of 
theoretically interesting models belonging to this class are: models with light leptoquarks; 
supersymmetric models with broken R-parity; TopColor models; some models based on 
non-Abelian horizontal gauge symmetries. 



A. Highly constrained models 

1. Minimal Supersymmetric Standard Model (MSSM) 

In a large part of the SUSY parameter space, the MSSM (for review and notations see, 
e.g., ||3^) is known to produce large effects on the radiative decay B ^ Xg'y, as well as on 
B Xsi~^ i~ (see [31,32|). The effects on B ^ XsVp have been studied in |32.3^. It was 



found that the contributions to the rate can be non-negligible only for tan (3 close to unity, 
while for increasing values of tan (3 the prediction rapidly converges to the SM value, regard- 
less of the values of the other SUSY parameters. Even for tan/5 ^ 1 it seems to be prob- 
lematic to enhance the rate up to the level observable at LEP, while keeping BR(i? X^ 7) 
within the experimental limits. This is due to the fact that the SUSY contributions to the 
657 vertex tend to dominate over those to the hsZ vertex. The SUSY corrections to the 
i? ^ 7 decay can be kept small while allowing for a large BR(i? z/), only for low 

values of tan /5, and by invoking large cancellations between the charged Higgs and chargino 
contributions. It is then conceivable that with a specific choice of several SUSY parameters, 
a "fine-tuned" MSSM could produce BR(i? ^ XgUu) close to ( p..4| ). However, we regard 
such a choice as unnatural. 
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2. Multi Higgs Doublet Models (MHDM) 



MHDM (for review and notation see, e.g.^ |31,3^) are severely constrained by B ^ -^s7, 
Z ^ bb, B ^ XcT u, and lepton universality in tau decays. The same is true for the 
more familiar two Higgs doublet models (2HDM), which represent a subclass of the gen- 
eral MHDM with natural flavor conservation. In a general MHDM the single parameter 
tan/3 of the 2HDM is replaced by three complex coupling constants, X, Y, and Z, which 
describe the Yukawa interactions of the lightest charged scalar with the down-type quarks, 
up-type quarks, and charged leptons, respectively. The new Z penguin diagrams present 
in these models are related to new photon penguins, and thus are severely constrained by 
B ^ Xg'y and cannot contribute significantly to B ^ Xsvv. A large enhancement not 
affecting i? ^ 7 and B ^ Xg fi~^ fi~ could arise for large enough values of the r z/7- 
Yukawa coupling via charged scalar box diagrams involving two external z/^-'s. These box 
diagrams are proportional to the combination Z^ (m^ + nis X^) ||3^. However, Y 
is constrained from Z ^ bb, the coupling Z from lepton universality in tau decays, and the 
product XZ from B ^ XcT u [p4tp!5| , p5| , implying together that also the contribution of the 
box diagrams to i? — > z/ z/ has to be very small. 



3. Left-Right Symmetric Models (LRSM) 

In these models (for review and notations see, e.g., |^) new heavy charged Wr gauge 
bosons coupled to right-handed currents are exchanged in a new set of electroweak penguin 
diagrams, which could give rise to deviations in the rate for B ^ XgV v . However, new 
photon penguins would enhance -B — > X^ 7 as well. The B —>■ Xg'j decay rate has been 
calculated in the minimal as well as in some non-minimal versions of the LRSM. Once the 
limits on Myy^ and on the Wl — Wr mixing angle are imposed (from direct Wr searches, 
Kl — Ks mass difference, etc.), the rate for i? ^ Xg 7 cannot differ from the SM prediction 
by more than about 50% The CLEO measurement of B Xs'~f further constrains 
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the parameter space of non-minimal models to regions where the SU{2)ji gauge coupling is 
small. In this region of parameters, the Z penguins are similarly suppressed, and therefore 
significant deviations of the B ^ XgV v decay rate from the SM prediction are excluded. 



B. Weakly constrained models 

1. Mixing of the b with new exotic Q = quarks 

It is well known that the presence of new heavy fermions with non-canonical SU{2) 
transformations (L-handed singlets and/or i?-handed doublets) mixed with the standard 
leptons and quarks would give rise to tree level FCNC in Z interactions |^|. In particular, 



the presence of new Q = isosinglet quarks, Dl and D^, as they appear for example in 
the 27 representation of Eq and in several superstring inspired extensions of the SM, would 
generate a hiSiZ vertex. New SU{2) doublets (^)_r ^^^^"^ with the i?-handed rf-type 

quarks would give rise to the new FCNC operator bRSnZ, which is absent in the SM. Both 
effects could appear simultaneously in the presence of a set of multiplets of mirror fermions 

It is easy to see why such flavor changing vertices are generated. For each L, R chirality 
state, the vector "^I ji = (^lo)^ ^ of the ordinary d° and new exotic D° quarks couples to the 
Z through the matrix of the isospin charges 

which, in contrast to the SM, are not proportional to the identity matrix. Then the isospin 
part of the neutral current for the mass eigenstates, '^l,r = Uj^ r^^l Ry which contains the 
flavor changing part of the coupling, reads 

>c„.ix = E [n 4 r n + ^R 4 r z, . (5.2) 

Here and henceforth the hermitian conjugate is understood. The strength of the flavor 
changing coupling kI ^ is given by 
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{1^3). (5.3) 



As a result, a mixing of the h with exotic rf-type quarks will induce a hqZ {q = d, s) vertex 
at tree level. We stress that a mixing is allowed only between states which carry the same 
quantum number of an unbroken gauge group. Therefore, for the corresponding photon 
coupling we have WQU oc I (with Q = — Namely, the electromagnetic couplings 
remain flavor diagonal. The same holds for the part of the Z coupling proportional to the 
electromagnetic generator. 

The effective interaction ( |2.1| ) arising from ( p.2| ) yields the coefficients 

^^'^=4cos2Vm|'''^'^ ^ ^'•^^i^§J''^'^''' ^^-^^ 
where the factor of 3 arises from the sum over the neutrino flavors. The current bounds 
on Ih'l r\ are obtained from the limit on B ^ /i+ /i~ ( [4.3|) (for q = s) and from — B^ 
mixing (for q = d) p9| , ^ . Numerically, they are comparable with the bound implied by the 



hmit O 

< 2.0 X 10-^ (5.5) 

In some mass matrix models, the flavor changing couplings generated by the ordinary- 
exotic fermion mixing are expected to be of the order of the ratio of the light to heavy 
mass scales (see however [^] for some interesting exceptions). We see that the limits 
on hsZ transitions are sensitive to h mixing with heavy particles up to a mass scale M ~ 
0{y/mf,ms/ k'"^), namely about 500 GeV. 



2. Models with large bsZ' vertex 

Quite in general, new fermions have to be present in models based on extended gauge 
groups (rank > 4), since they are needed to ensure the absence of anomalies in the new 
gauge currents. Eq models are a well known example where new SU{2) singlets Dl and Dr 
are present, together with new neutral gauge bosons corresponding to the two additional 

20 



Cartan generators of the group. We generically denote the new gauge bosons as Z\ which 
are coupled to the fermions through charges Q' . A mixing with new ordinary quarks {i.e., 
with conventional SU{2) quantum number assignments) would not affect the couplings of 
the Z. However, as the matrix of the Q' charges is in general not proportional to the identity, 
H = (^I Q'l r Ua)ij^j "will not vanish, thus inducing a bsZ' vertex. 

In spite of the 1/M|, suppression, the Z' mediated FCNC are expected to be as large 
as the corresponding transitions induced by Z exchange This is because while in 



general the flavor changing mixings affecting the Z couplings are suppressed as the ratio 
of the light to heavy fermion masses, no analogous suppression is expected for the mixings 
between fermions of the same isospin, which affect only the Z' couplings. The absence of the 
suppression in the mixing can compensate for the M'^/M'^, suppression of the Z' relative to 
the Z amplitude, implying that the coefficients describing the Z and the Z' flavor changing 



effective interactions can be comparable in size . 

The analysis of the NP effects on an enhanced bsZ' vertex parallels quite closely that 
of the large bsZ models. However, in the former case the enhancement of the B ^ XgU u 
decay mode relative to B XgC,^ [~ is not as large as in the latter one. Since and ^L 
appear in an SU{2) doublet, they have the same charge. Then the ratio of the rates of 
the two processes is given by 



(5.6) 



In most models, the Z' couples universally to all three generations, and then the above 
ratio cannot be larger than 3. (This upper bound corresponds to the case when the right- 
handed leptons are almost decoupled from the Z' , namely Q/j ~ 0. This can happen 
in Eq models in which the Z' arises as a particular combination of the two additional 
f/(l)' generators.) As a result, these models are significantly more constrained by the limit 
on B —>■ Xs fJ^^ fi~ ( [4.2| ) than the models with a large bsZ vertex. The UAl limit implies 
BR(i? ^ XgUu) ^ 1.8 X 10^^, which is only marginally within the reach of the expected 
LEP sensitivity (|r|). 
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The expressions for the Cl,r coefficients are similar to 



Cl,r - ^l,r{Q') =^ C^^R - 3 



m: 



2 



Mi, 



bq 1 2 



q=d,s 



(5.7) 



Here i^'lr describe the strengths of the FCNC couphngs of the Z' (the analogs of ( ^.3|) ), and 
we defined 

^L,i^(Q') = Q\v) [Q\dL,R) - Q\DL,n)] , ^9 = ^ • (5-8) 

For any particular Z' model, J^l,r and are known. Then the limit ( |1.2| ) implies constraints 
on the ratios K^-'^/Mf,, typically about a factor two weaker than the order 10~^ limits from 
the present bound on i? ^ /i^ 

Thus, in this class of models the largest allowed effects can hardly enhance the 
B ^ XsV v decay rate up to the estimated LEP sensitivity (|1.4| ). However, the sensitivity 
of the B XgV v decay rate to Z' effects can be larger in a class of unconventional Eq 



models In these models, the different generations are embedded in three fundamental 
27 representations in a generation dependent way, implying different charges for the left- 
handed leptons of the different families. In such a scenario the muon can be weakly coupled 
to the Z' without implying the same for and i^e, and then the ratio of the B ^ XgU u 
and B X^yU+yU^ decay rates can exceed the previously derived limit of 3. For the uncon- 



ventional Eq models of Ref. |^ we find that the ratio in ( [5. 6] ) can be as large as 5. In this 



-4 



case the constraint from B ^ Xs fj,^ fj, ( [4.2|) still allows BR(i? — > z/) up to 2.9 x 10 

3. Large tcZ effective coupling 

In several NP models, new sources of FCNC are naturally suppressed, as they are related 
to ratios between the masses of the fermions involved in the flavor changing transitions and 
some large mass scale. Due to the large value of the top mass, such a suppression might not 
be effective for flavor changing transitions involving the top quark ||4^,^ , and theoretically 
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the presence of a large tcZ vertex is indeed an open possibility (see [^6|-|48[], and references 
therein). This can be the case in models which predict new dynamical interactions of the 
top quark p9| . |50[] , in MHDM without natural flavor conservation or in the presence of 



mixing with new Q = ^ isosinglet quarks P8| , p2| , p3 

To date, the couplings of the top quark have not been directly measured, and a large 
tcZ vertex is not constrained by the measurement oi B ^ Xg 7, neither by — B^ mixing. 
However, an anomalous tcZ coupling will give new contributions to the effective hsZ vertex 
46| , |48|| . Contributions to the B XsV v decay would arise from diagrams involving loops of 



charged bosons and of unphysical 0^ scalars, with an insertion of the tcZ coupling on the 



fermion line Since these diagrams are not CKM suppressed (they are proportional 

to Vj^ Vcs), there is no additional suppression beyond the loop factor. 

The case of a tree level tcZ vertex induced by a mixing with new <5 = | isosinglets was 



analyzed in In this particular case the vertex arises at tree level. Therefore, after the 



underlying theory is fully specifled, the computation of the new penguin diagrams can be 
performed in detail, yielding a flnite result |^S|. It was found that the new contributions to 



the hsZ effective vertex are always smaller in absolute value than the SM contributions, and 
opposite in sign. As a result, they interfere destructively thus lowering the expected rates 
for B ^ Xsd-^ and B XsV v . Thus, an anomalous tcZ vertex induced by mixing with 
Q = I isosinglets cannot be constrained by an upper bound on the branching ratio. 

If the tcZ vertex is an effective one, and the underlying theory is not specifled, the 
expression for the loop induced hsZ vertex by itself is formally divergent, and it has to be 



regulated. We take the result of the computation of the flavor changing penguins from |^ 
and we substitute the function resulting from the flnite loop integration of that speciflc case, 
with a regulator log(A^/m^). We further assume that new effects dominate over the SM 
contributions. This yields for the induced hqZ vertex 



cose 



w 



, . 2n (v-;,4^K,)iog— 



(5.9) 



where parameterizes the strength of the tcZ coupling. For the dimensionless coefficient 
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defined in ( p.lO| ) we get 



a 



4 V47rsin^6'v^ 



Kt\ 1o£ 



q=d,s 



eg I 



(5.10) 



Assuming A ~ 1 TeV, \Vtb\ ~ \Vcs\ ~ 1, \Vcd\ ^ 0.22, we get from 



log^ < 0.55 



141 ^0.16^ 



(5.11) 



wliicli is comparable witli tlie limit given in once the differences in the normalization 
and in the overall coefficient in (5.91) are accounted for. 



We have chosen to present our bound in the form (|5.11| ) for uniformity of notations, 



and to allow for comparison with However, not to mislead the reader, a remark is in 
order, parameterizes an effective vertex, and thus if the heavy physics decouples from low 
energy, this parameter must be proportional to inverse powers of the NP scale A. Therefore, 
even if this is not apparent from the notation we used in ( |5.9|) and from the one used in 
P6| , in the K ^ oo limit the hsZ coupling approaches its SM value, as it should, in order 
to account for the decoupling of heavy physics. 



4.- Fourth Generation 



Extensions of the SM including a fourth generation are still an open possibility 



and there are theoretically motivated models in which the additional neutrino is naturally 
heavy |5^, thus escaping the LEP limit on the number of light neutrino species. It has also 



been shown that SUSY models with four generations are consistent with unification, and 
by imposing this requirement some conditions on the t' , 6', r', and v' spectrum have been 
derived [^]. In these models large CKM mixings with a heavy t' may enhance the Z and 



photon penguins 

The constraints implied by the CLEO measurement of 5 — 7 were presented in 



After constraining the new mixing angles by the limits on deviations of the 3x3 CKM 
matrix from unitarity, it was found that a new t' quark in the mass range 200 — 400 GeV is 
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still consistent with experimental data. The possibility of having a measurable enhancement 
of the Z mediated FCNC decays, while keeping B Xg^ close to the SM prediction, relies 
on the fact that the effective 657 vertex has a logarithmic dependence on mt/, while for hsZ 
this dependence is quadratic. The dependence of the B Xgvv rate on mf is numerically 
similar to that of -B — > Xgfi^fi^. The different contributions of the box diagrams tend 
to slightly enhance the neutrino decay over the charged lepton mode. We find that the 
constraint ('14.2D still allows 



BR{B ^ X,vv) < 3.7 X 10"^ . 



(5.12) 



In the presence of a heavy t' quark, the SM coefficient ( p.5| ) is modified according to 

2 



1 + 



(5.13) 



ViVts Xoixt) 

where, for simplicity, we assumed \Vt's\ ^ and \Vts\ ^ l^tdl- From the limit on 

B ^ XsV v ( |1.2D we get the bound 



which is numerically similar to the limit implied hy B ^ Xsd,^ 



(5.14) 



5. Anomalous WW Z couplings 

The SU{2)l x f/(l)y gauge symmetry of the SM fixes the dimension-4 operators that 
describe vector-boson self-couplings. This symmetry should be respected by any low en- 
ergy effective theory, independently of possible NP at energies above the electroweak scale. 
NP may still signal itself in low energy experiments through higher dimensional operators, 
suppressed by inverse powers of the NP scale A. While dimension-6 operators modify the 
WW'-f and WW Z vertices identically, this is no longer true for dimension-8, and higher 
dimensional operators [^. Thus, as long as dimension-6 operators dominate the possible 
deviations from the SM, the experimental measurement of -B ^ 7 implies that also the 
hsZ coupling must be close to its SM value. 
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In general, dimension-6 operators, which are suppressed by two inverse powers of the NP 
scale, are expected to dominate over dimension-8 operators, suppressed by A~^. However, 
in any perturbative underlying theory the dimension-6 operators can only arise from loop 
diagrams, and therefore they have an additional loop suppression factor of l/167r^. In 
contrast, dimension-8 operators can also arise at tree level [BUI. Thus, if the scale of new 



physics is below about 2 TeV, then dimension-8 operators can dominate over the dimension- 
6 ones. In this case, the measurement of i? ^ 7 implies no direct constraints on the bsZ 
vertex. 

The effects of anomalous WW Z couplings on rare K and B decays were studied in 



3| . In general, seven coupling constants parameterize the WW Z interaction ||6J] . In low 
energy processes, when the external momenta can be neglected, only gf and contribute 
T|,|62[|. A recent CDF measurement implies that Agf = gf — 1 has to be small. 



Therefore, for simplicity we assume gf = 1, and we study only the effects of g^ . The Z 
penguin ( |2.6| ) contribution is modified according to pT|-|63|]p| 

Xoix,) Xoix,) + gl cosl Wix,) , W{x) = ^"^^^ "^'^J"^^,^'''^^ • (5.15) 



The strongest published bound on g^ is obtained from the measurement of Kl jj,^ pT 



As we have discussed, because of the large long-distance contribution to this process, such a 
bound is not very reliable. We prefer to quote the bound resulting from BR(i? Xg yU+ 
For nit = 180 GeV we found that the limit on this decay mode implies —8.3 < g^ < 3.7. 
This constraint allows for BR(i? — > z/i/) up to the value in Eq. (^4.3| ). Our bound ( |1.2|) 
implies 

-8.6<^f<4.1. (5.16) 



We use the sign convention of |63| for , which differs from that of |61,62|. 
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6. Extended Technicolor (ETC) 



FCNC processes impose very strong constraints on technicolor models. Either "tradi- 
tional" ETC models with a minimal set of interactions necessary for third family quark mass 
generation, or ETC models which incorporate a techni-GIM mechanism, yield a. B Xs'j 
rate at most slightly larger than the SM rate While at leading order the bs'j coupling 
is not affected, these classes of ETC models typically induce large flavor changing Z boson 
couplings of the form si'^^hi . Techni-GIM models yield an even larger enhancement of 
4-fermion interactions of the form Ol in ( |2.2| ). 

It was subsequently noted that ETC models with a techni-GIM mechanism also 
predict the B ^ Xs fi^ fi^ decay rate about a factor of 30 above the SM, violating the 
experimental bound ( [4.2|) , unless significant cancellations occur between various contribu- 
tions. However, "traditional" ETC models predict only about a factor of 4 enhancement 
of B —* Xg i~ over the SM. Therefore, these models are not yet excluded, and yield a 
similar enhancement for the B ^ XgU u decay rate, which could be within the reach of the 
expected LEP sensitivity (p^). 

C. Unconstrained models 

1. Light leptoquarks 

Leptoquarks (LQ) couple directly leptons to quarks. Such particles appear in several 
extensions of the SM. A comprehensive analysis of the experimental constraints on the 
LQ couplings has been given in |68|, and is summarized in Table 15 of this reference. After 



summing over all the possible neutrino flavors in the final state, it turns out that the existing 
limits on the LQ Yukawa couplings allow for BR(i? ^ Xquu) up to (9(10%). For the 
relevant couplings, our new limit (|1.2|) imposes much stronger constraints than the existing 
ones. Several types of LQ are possible, and we adopt here the notations of [^. LQ can 



be scalar (S) or vector (V) particles, and can belong to different SU{2)l representations. 
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The SU{2)l singlets, doublets, and triplets are labeled with the lower index 0, 1/2, and 1, 
respectively. 

The following LQ can mediate the B ^ XsV v decay: 

^0, ^1/2, ^1, ^1/2, ^1. (5.17) 
Written in components, the relevant scalar and vector terms in the interaction Lagrangian 



are 
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C-LQ - - Afg ql + Xiq'^ Qr v\ 5*1/2 - Ql ^1 

+ a5^^ q'kl, vl Vr/l + Ag QLl, ylVr\ (5.18) 

where q = d,s,b, and it is understood that only the charge | component of the SU{2)l 
multiplets appears in (|5.18|) . For simplicity, we assume that all the \iq couplings are real. 



Integrating out the LQ fields and Fiertz transforming, yields for the coefficients of the 



effective four-fermion interaction ( p.l|) induced by ( 5.1 



Ct^R = VLQ^^. {q = d,s), (5.19) 

where rjiQ = 1/2 (1) for scalar (vector) LQ. The coupling Cl is generated through Sq, Si 
and Vi exchange, while Cr appears from S'l/2 and V1/2. As different types of LQ can exist, 
both Cl and Cr can be simultaneously present. If in addition the LQ carry some generation 
index, cancellations between different generations of LQ are also possible, and then the limit 
( |1 . 2|) constrains only the total LQ-mediated rate. For simplicity, we will restrict ourselves 



to the case when only one type of LQ is present, and it does not carry any generation index. 
For scalar LQ the limit ( |1.2|) implies the following new bounds 

while for vector LQ 
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These bounds are much stronger than the existing hmits . 



Other LQ couphngs involving the hght fermions of the first and second generations are 
constrained by the existing experimental data to be much smaller than the limits ( |5.20|) 



and ( |5.21| ) ||69|| . However, if we have to learn a lesson from the hierarchy in the fermion 



Yukawa couplings, then it seems natural to expect a large hierarchy in the LQ couplings to 
the different generations as well. This is the case, for example, in models that explain the 
quark and lepton mass hierarchies as originating from horizontal symmetries. If LQ exist, 
in these models they couple more strongly to the third generation fermions Since three 



third generation fields and only one from the second generation participate in the process 
B ^ X s Ut- , any improvement in the search for the B ^ XgV v decay would represent an 
important test of these models. 

2. SUSY with broken R-parity 

In SUSY models it is usually assumed that R-parity is a good symmetry. However, this 
is not necessarily the case, and one can construct SUSY models with broken R-parity. We 



concentrate on the MSSM without R-parity ||7l|. Extra trilinear terms are allowed in the 



superpotential, and some of them can give rise to a large enhancement of the B Xg v v 
decay rate. Denoting by L^, Q\ and d\ the chiral superfields containing respectively the 
left-handed lepton and quark doublets, and the right-handed down-type quark singlets of 
the i-th generation, these terms read 

Ty^ = A^,LlQi4, (5.22) 

where, for simplicity, we assume the A^^^ couplings to be real. Omitting terms involving the 
Ml and fermions which are not relevant for the tree level h —>■ qviVj transition, the Yukawa 
interactions of the R-parity breaking Lagrangian generated by ( |5.22| ) are 



= K,k [ di (d'n ^I) + 4* di) ] . (5.23) 
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The exchange of da and di squarks gives rise to the effective four-fermion interaction 
l]) responsible for B XgUii/j, with the coefficients 

^f = E%^^ ^f=E%$^, i9 = d,s). (5.24) 
k ^'"d*^ k ^'"jfc 

We note that, in contrast to the LQ case, here both the Ol and 0/j operators are necessarily 

present, due to the simultaneous appearance of the dL and dR scalar superpartners of the di^R 

quarks. For simplicity, we neglect possible cancellations among the different # {k = 1, 2, 3) 

exchange amplitudes. Then the bound ( |1.2| ) implies the following limits on the product of 

R-parity violating couplings 

AUA5.<i.ixio-'(^)\ a;„a;,3<i.ixio-(^)\ (5.25) 

These limits represent the strongest constraints on the product of couplings involving 
the third generation, such as A3g3A333, A33gA333. As in the LQ models, these couplings are 
expected to be particularly large in models that relate their size (relative to the couplings 
involving the lighter generations) to the fermion mass hierarchy 173 . 



3. Top Color models 



TopColor models attempt to explain the large value of the top mass through the 
dynamical formation of a ti condensate. In these models, the basic assumption is that new 
dynamics strong enough to form chiral condensates is effective for the third generation, 
which therefore is treated differently from the first two. As a consequence, TopColor models 
have peculiar implications for the phenomenology involving the third generation, such as 
top and bottom production []73[ and the Z ^ bb decay rate [|7S|. In particular, they are 
expected to yield large effects in FCNC B decays involving third generation leptons |7B| , |77| . 

Several models can be constructed along these lines, and we concentrate on the one 



studied in In this model the gauge symmetry breaking structure is 



SU{3)i X [/(l)i X SU{3)2 X f/(l)2 ^ 5f/(3)QCD x f/(l) 



Y 



(5.26) 
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Here SU{3)i x f/(l)i couples only to the third generation while SU{3)2 x f/(l)2 couples 
only to the first and second generations. The quantum numbers under these groups coincide 
with those under the usual SU{3)qcd x U{1)y- The breaking into the SM group ( |5.26| ) 
is induced by a (tt) condensate, which is generated at the 1 TeV scale when the SU{3)i 
coupling becomes strong. The initial symmetry is larger than the electroweak gauge group, 
and this implies the existence of new massive gauge bosons corresponding to the additional 
broken generators: a color octet (topgluons) and a singlet Z'^. We concentrate on the 
Z' boson since it can mediate B ^ XsV v decays. The couplings of Z' to the fermions are 
given by IT^l 



(5.27) 

with Ql = {u,d)L, Ll = {i,h')L, gi ~ 0.35 is the f/(l)y coupling constant, and z = 1,2,3 is 
a generation index. 

The most important difference between this model and the usual Z' models arises from 
the fi{d') factor, which enhances the strength of the third generation couplings with respect 
to the first and second generations. One has /i,2(^') = tan^^' for the first two generations 
and fsiO') = —cot 6' for the third generation. In general cot^' ^ 1 is expected, in order 
to ensure that the condensate forms in the top direction After integrating out the 



Z' boson and rotating the dL,R quarks into the mass basis by the unitary matrices Ul,r, 
the interaction Lagrangian ( [5.27] ) gives rise to the following coefficients for the effective 
four-fermion interaction (12? 



r^qn _ 1 al ,bqr _^ ~2 _ /sill^ M| > 



12 M|, ' V 6 Ml 



2 



1=1,2,3 
q=d,s 

2 



q=d,s 

Here k""' = J2j{U^j Ujq)fj ~ —{U^sUsq) cot 9' gives the flavor changing Z' couplings to the 
quarks (the L and R chirality labels for k'^^^ and for Ul^r are understood). 
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A too large value of cot 9' would lead to a spontaneous breaking of the chiral symmetry 
for the tau lepton. This implies the bound |177| , |78 



coi^e'<^ =^ |cote'|< 25. (5.29) 



9l 

Inserting the values of the various hypercharges given in ( p.27| ) into Eq. ( p.6| ), we find 

BR{B ^Xsvv) ^ n + /I + /I ^ cot^ 
BR(5^X,/i+/i-) 5/1 ^ 5 ' ^ ' 

where we assumed that NP effects dominate the decay rates, and thus we neglected the 



SM contribution. As is apparent from ( |5.3CI| ), due to their different dependence on 6*', a 
measurement of both B ^ Xg fj^^ and B —>■ Xsiy jy would allow us to separately determine 
the value of this parameter. The current limit on S — » z/ in this model is obtained by 
combining Eqs. ( [5.29| ) and ( |5.30| ) with the experimental upper bound on B ^ Xs fi^ ^~ 
( [4.31 ). The B ^ XgV v rate is constrained rather weakly by this process, which still allows 
for a branching ratio up to about 10%. The bound ( |1.2| ) gives the new stringent limits 

Icot 9' ^'!"^\ < 5.4 ( ' , Icot 6' z.;^"! < 2.7 



, Icot 9' K^'^l < 2.7 — ^ . (5.31) 
ITeV/ ' ' VlTeV/ ^ ' 
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4- Horizontal gauge symmetries 

Attempts to explain the hierarchical pattern of fermion masses and mixings by some 
underlying dynamical interaction, are often based on broken horizontal gauge symmetries 
[|79| , |80|j55| . In the non-Abelian case, the fermions of different generations are assigned to 
some irreducible representation of the horizontal gauge group. Flavor changing transitions 
occur, suppressed by the masses of the heavy horizontal gauge bosons. 

The fermion mass pattern is generated through the so-called universal see-saw mechanism 

-—1/2 

|81|, in which the fermion masses are suppressed from their natural scale Gp , to the 
observed values by inverse powers of some large mass scale, with no need to fine tune the 
Yukawa couplings [ [75| , |5n| . Various mass scales are associated with different stages of the 



horizontal gauge symmetry breaking in such a way that the heavier fermions of the third 
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generation couple to the lightest horizontal gauge bosons. As a consequence, rare FCNC 
transitions of the third family are "naturally" enhanced. 

To illustrate the main features of these models, we consider as a simple example a model 
based on the horizontal gauge symmetry group SU{3)h x U{1)h [0- The standard q 
quarks are assigned to the (3, 1)h representation of SU{3)h x U{1)h, while new isosinglet 
Q quarks are in (3, —1)h- The breaking of the horizontal symmetry is achieved through the 
non- vanishing VEVs C,a, of a set of SM singlet horizontal scalars belonging to the (3,2)// 
or (6, 2)h, which also give large masses to the heavy isosinglet quarks. The mass terms for 
standard quarks are provided by the VEV of the SM Higgs, v = {2^/2GF)~^^'^ — 175 GeV, 
and by an additional VEV rj of a real scalar isosinglet belonging to (1, 0)h- Then, for the 
fermion mass terms we have 

Cm--QIM {0^J Qi + (fLQ'RV + Qi <fR V , (5.32) 

where i,j = 1,2,3 are flavor indices. We note that due to the assumed horizontal gauge 
symmetry, the second and third terms in ( p.32| ) involving the light fermions are the same 
for the three generations. 

By assumption t>, r/ <^ , and thus the quark masses are generated via a see-saw like 
mechanism, yielding for the mass matrices rriq ~ vrj /iA{^)~^. All information on the quark 
mass hierarchies and mixings are contained in and depend on the hierarchical struc- 

ture of the VEVs ^Q. It is assumed that ^ ^2 ^ ^3 in order to reproduce the generation 
hierarchy pattern. For example, together with additional smaller induced VEVs contributing 



to A4, they can be chosen to form a Fritzsch structure yielding 



6 : 6 : ei ~ 1 : J— : J— ^ ~ 1 : 190 : 3300 . (5.33) 

y m„ \ ruu rriu 

The horizontal VEVs hierarchy induces the breaking chain 

5f/(3)H X f/(l)H ^ 5f/(2)H X f/(l)'^ ^ f/(l)'^ ^ /. (5.34) 

In the first stage, four flavor changing gauge bosons carrying family "charge" plus a com- 
bination of the three neutral generators acquire a large mass of order ^1, leaving a residual 
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SU{2)h X (/{lYfj symmetry unbroken, which acts only on the second and third generations. 
The four "charged" bosons couple the first family fermions directly (namely, not through 
small mixing angles) to the fermions of the second and third generations, and in particular 
they induce AS = 2 effective operators. The requirement that these operators will not 
generate unacceptably large contributions to the — mass difference implies for the 



scale of the first breaking 



^1 > 3 X 10^ TeV . 



(5.35) 



In the second stage, two charged bosons and a second neutral combination acquire a mass 
of order ^2- These bosons give direct contributions to the transition B — > XsVtI^hi and 
as a consequence our limit (|1.2|) implies a lower bound on ^2 (see Eq. ( ^.38|) below). The 
remaining local symmetry U{1)"^ that acts only on the third generation is broken at a scale 
^3 much smaller than ^1 and ^2- For example, if we assume the hierarchy ( |5.33| ), the bound 
( p.35| ) implies the rather weak constraint ,^3 > 0.9 TeV. 

We see that this scenario has the natural implication of a new X° boson corresponding 
to U{iy^ in ( |5.34| ), coupled only to the third generation, while due to the hierarchy (^ 



the effects of the heavier horizontal bosons are much more suppressed. Due to the mixing 
between the quark mass eigenstates, X° will also mediate flavor changing transitions. The 
neutral current interaction Lagrangian describing these transitions is analogous to (|5.2| ), and 
the terms relevant to the process we are interested in are 



bq- 



1^0 



(5.36) 



where kl^ji = Ul,r^P3 Ul,r with P3 =diag(0, 0, 1), and the matrices Ul^r rotate the dL,R-type 
quarks mass eigenstates. Eq. (|5.36|) yields the coefficients 



<-'L,/? — 72 '^L,R 
S3 



C 



L,R 



- y 

^3 q=d,s 



bq 1 2 

l,r\ 



(5.37) 



If the FCNC mixings k^^^ are very small, then in spite of the larger mass suppression 
the direct transition B — > z/^ z/^ mediated by the "charged" SU{2)h bosons of mass 
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Mx± = 5'_H"^2 (see ( p.34| )) can give contributions to i? — > X^z/z/ which are competitive to 
those mediated by The rate for the mediated decay B ^ XgU^ v^, becomes as large 
as the rate for B z/^ z/^ mediated by X° when r\ ~ ill ^2- Due to the different 

final states, the two amplitudes do not interfere, and the coefficients Cl,r for the second 
case follow easily from ( |5.37| ) by the substitution i^^l jiHl — > l/^l- ^™ limit ( |1.2| ) implies 
the two bounds 

A different scenario which also predicts a new U{1) bosons coupled only to the third 
family fermions was presented in Similar to TopColor models, in this scenario the 

large value of the top mass is explained by a dynamical scheme which implies a large isospin 
breaking for the t and h masses. However, that does not feed back directly into the W and Z 
masses. The U (1) gauge boson coupled to the third generation is the remnant of the breaking 
of some large non-Abelian semi-simple gauge group which embeds the HyperColor interac- 
tion responsible for the formation of the isospin breaking condensate. The phenomenology 
of such a gauge boson, including the consequences on LEP physics due to its mixing with the 
Z, have been extensively studied |8^. Recently it was also suggested that the B Xgiy u 
decay rate could be largely enhanced because of the new FCNC contributions, and could 
even approach the rate of semileptonic B decay |Q . Such a mechanism for increasing the B 
width was proposed as a possible solution to the claimed discrepancy between the observed 
experimental value of the semileptonic B branching ratio and the theoretical predictions. 
The limit ( |1.2|) obviously rules out this possibility, as well as any other mechanism attempt- 
ing to achieve a sizable increase of the B width through an enhancement of decay modes 
associated with large missing energy. 

The interaction Lagrangian for the X^ boson of this model, written in the fermion mass 
basis, is identical to the Lagrangian in Eq. ( [5.361 ). However, there is now one additional 



constraint for the X^ mass ||8^ 
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Using Mx = Qh^^i it is then straightforward to derive from (|5.37| ) the expression for the 
C L,R coefficients 

cIr = Ii E \<R?- (5-40) 

q=d,s 

The NP parameters k^j^ are denoted by and — (z = 1, 2) in Ref. [Q. It was speculated 



that P + P ~ 30 |Vcbp could account for the possible discrepancy in the B semileptonic 
branching ratio [Q. Our limit ( [L.2D implies a bound two orders of magnitude smaller, 
I'^L P + l^flP < 2 X 10""^. For q = d, the limit \k,''l — k,r\ < 2 x 10~^ can be derived from 
Bd — Bfi mixing [Q. However, there is no similarly strong limit on the individual 1^5^^! 
parameters. From (|1.2| ) we obtain 

< 1.4 X 10-2 (g = d,s). (5.41) 



VI. SUMMARY AND CONCLUSIONS 

In this paper we discussed the derivation of the first bound on the inclusive B ^ Xg u u 
decay rate, using the large missing energy tag in B decays at LEP. We studied in detail 
the theoretical ingredients needed to carry out such an analysis, and we found that the 
overall theoretical uncertainty is small. Therefore, this decay mode is well suited to search 
for physics beyond the SM. We translated the ALEPH bound on the B ^ tv branching 
ratio, which resulted from a search for B decays with large missing energy, into a limit on 
the B ^ XgU u branching ratio. To derive a numerical limit, we had to make a number of 
conservative and simplifying assumptions. Thus, the resulting bound is weaker than what 
a dedicated experimental analysis will be able to achieve. Our conservative upper bound is 

BR(5 ^ X,z/z/) < 3.9 X 10-\ (6.1) 

which is less than one order of magnitude above the SM prediction. We estimated that 
using the full LEP-I data sample, the LEP collaborations may be able to set a limit of order 
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(1 — 2) X 10~^. Due to the theoretical interest in the B XsV v decay mode, we think it 
is important that the LEP collaborations will perform a dedicated analysis of this process. 

We studied a variety of new physics models. After discussing the constraints from existing 
experimental data, we divided the NP models into three classes, according to increasing 
allowed values of the the B ^ X^v v branching ratio. 

To class (A) belong those models in which the B ^ XsV v branching ratio is already 
constrained to be below the expected sensitivity of the LEP experiments. 

Class (B) contains the models which allow for an enhancement of BR(S ^ z/ z/) up to 
values that will be observable at LEP. These models are listed in Table |I[ together with the 
maximal B XgV v branching ratio they allow for, after the existing constraints are taken 
into account. These models naturally evade the constraints imposed by 5 — * Xg'-y. The 
limits on the relevant NP parameters which we obtain from the bound on the B ^ XgU u 
decay rate, are numerically close to the limits provided by the bounds on the inclusive and 
exclusive B ^ Xsi~^ i~ decays. However, our bounds are more reliable, as the B —>■ Xg u u 
decay is theoretically cleaner. We expect that in the near future, new results emerging from 
CLEO, CDF and (hopefully) from LEP can compete in further constraining these models. 

The most interesting models for our investigation belong to class (C). For the models in 
this class, the bound on the B ^ XsV v branching ratio implies new constraints that are not 
matched by other existing experimental data. A generic feature of these models is that they 
yield a natural enhancement of the FCNC processes involving third generation fermions, 
without conflicting with other constraints. We derived new limits on the couplings of the 
third generation fermions in models with leptoquarks, and in supersymmetric models with 
broken R-parity. Our bounds also imply stringent constraints on models in which new gauge 
bosons are coupled dominantly to the third generation, such as TopColor models, models 
based on non-Abelian horizontal gauge symmetries, and other models attempting to explain 
dynamically the large value of the top mass. The new bounds on the parameters of these 
models, implied by the limit on B ^ Xg u i?, are summarized in Table |T[ 

In future B factories much larger data samples will become available. While - to our 
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knowledge - no detailed study concerning the possibility of measuring the B ^ XgV V decay 
rate at B factories has been carried out, we hope that it will be possible to perform such a 
search. A precise measurement of the B ^ Xsfi' decay rate would provide a very reliable 
means of directly determining the iVt^l element of the CKM matrix. This would allow for 
an important new test of the unitarity of the CKM matrix. If deviations from the SM 
predictions will be detected in rare B decays, it will also be important to measure as many 
decay modes as possible. The pattern of deviations from the SM predictions in different 
decay rates would help us to distinguish between various possible NP scenarios. 
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TABLES 



Model 


Allowed branching ratio 


SM 


0.5 X 10-4 


FCNC Z 


3.5 X 10-4 


FCNC Z' 


1.8 X 10-4 


Unconventional 


2.9 X 10-4 


Anomalous tcZ vertex 


3.5 X 10-4 


Fourth generation 


3.7 X 10-4 


Anomalous WWZ 


3.5 X 10-4 


ETC 


3.5 X 10-4 



TABLE I. Summary of the models belonging to class (B), which allow an enhancement of 
the B Xg V V decay rate up to a level observable at LEP. The standard model (SM) prediction 
for BR(5 Xgvv) is given for comparison as the first entry. For each model listed in the first 
column, the second column gives the maximal B ^ Xsvv branching ratio allowed by the existing 
experimental data. 
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Model 


IVpw hoi in (is 


LQ: 5*0, 5*1/2) 


A.A,3< 1.1 X 10-3 (^-g^ 




A,,A,3<5.7xlO-^f ""^Q y 
^ VlOOGeVy 


SUSY without R-parity 


A:fc,A;.fc3 < 1.1 X 10-3 ^ 
^ V 100 GeV / 


Top Color 




Horizontal gauge symmetry 


|4«^|<6X10-^ (-rS^)^ 
V 100 GeV/ 

6 > 6TeV 


Model of Ref. [S3,3| 


< 1.4 X 10-2 



TABLE II. Summary of the models belonging to class (C), for which essentially no constraints 
on the allowed B Xg v v branching fraction existed to date. For each model listed in the first 
column, the second column gives the bounds on the relevant model parameters implied by the 
new limit BR(i?^XsZ^P) < 3.9 x 10"^. The indices A: = e,/i,r and q = d,s correspond 
respectively to the neutrino and quark flavors in the final state. 
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